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Abstract —In this paper we have obtained approximate 
solutions of a wave equation using previously studied 
method namely perturbation-iteration algorithm (PIA). 
The results are compared with the first and second order 
difference scheme solutions by absolute error. In 
addition, to prove the effectiveness of the method, we have 
presented some graphics and tables. 

Keywords—Initial value problems, wave equation, 
perturbation-iteration algorithm, difference schemes, 
error analysis. 


I. INTRODUCTION 

Partial differential equations have been used to 
explain many phenomena in different science and 
engineering branches such as mathematical biology, 
physics, image processing, quantum mechanics, fluid 
flow, viscoelasticity and so on. Therefore, to understand 
and explain physical interpretation of the problems arise 
in the above-mentioned fields, a considerable effort has 
been achieved and numerous methods have been 
proposed to obtain both numerical and analytical 
solutions of the partial differential equations. These 
methods include, Adonian decomposition method 
(ADM)[l-2], variational iteration method (VIM)[3-4], 
homotopy analysis method (HAM)[5-6], homotopy 
perturbation method (HPM)[7-8], finite difference 
method (FDM)[9-11], differential transform method 
(DTM)[12-13], etc. 

Alongside these methods, a perturbation-iteration 
method, namely perturbation-iteration algorithm (PIA) 
has been proposed by Aksoy and Pakdemirli in 2010 [14- 
15]. In this paper we implementPIA to obtain some 
approximate solutions of a wave partial differential 
equation with initial conditions. Obtained results are 
compared with the known exact solutions and the 
solutions obtained by the finite difference method via first 
and second order difference schemes. The findings are 
satisfactory and the present method produces highly 
approximate results even for a few iterations. 

II. BASIC IDEA OF PIA 

In this section we introduce some fundamental points of 
the PIA. 

Take the wave partial differential equation: 

F (u tt ,u t ,uu xx ,u,E) = 0 (1) 


where u = u(x, t) and £ is a small perturbation parameter 
that will be inserted to the equation later. The perturbation 
expansion with only one correction term is 
tt n+ i — u n + £ (u c ) n (2) 

Replacing Eq.(2) into Eq.(l) and writing in the Taylor 
series expansion with first order derivatives only gives 
F((u n ) tt , (u n ) xx , u n ,0) 

+ i'u 1 1 (.u n \, (u n ) xx , u n , 0) £ ((it c ) tt ) n 

+ F Ut Cu n ) xx , u n , 0)£((u c ) t ) n 

4” Pu xx ^Tt-n)tt’ (tinX' )xx’ u n> 0)£((li c ) xx ) n 

+ F u (u n ) t , (u n ) xx , U n , 0)£(li c ) n 

T F e C , (tt n ) xx i tt n , 0)£ 0 (3) 


or shortly, 

F 

J + ((.u c ) tt ) n F Utt + ((u c ) t ) n F Ut 


where 
and F Un 


+ (( u c ) xx ) n F Uxx+ (.u c ) n F u +F e = 0 (4) 

_dF _ 0F _ dF _ dF 

Utt du tt ‘ Utt du tt ‘ Utt du tt ‘ Utt flu tt 
dF 

3u t t 


hi the expansion, allof the derivatives are calculated at 
£ = 0. Opening with the initial assumption u Q (x,t ), in 
the first step (u c ) 0 (x, t) is determined from Eq.(3) and 
then subrogated into Eq.(2) to obtainii 0 (x, t).The iteration 
procedure continues until a desired solution is obtained. 


III. NUMERICAL RESULTS 

Consider the following wave equation [16] 

/6t 4 + 4t 2 — 2 x + t \ 

u tt ~ (* + t)u xx = ^ (l+fc 2)4 + ) sinW ^ 5 ) 

for 0 < t < 1, 0 < x < n given with the initial and 
boundary conditions 

u(x, 0) = sinGc) ,it'(x, 0) = 0, 0 <x<n 
u(0, t) = u(n, t) = 0,0 < t < 1. 

The known exact solution of the problem is 
1 

u(x,t ) =-rsinGc). (6) 

1 + t 2 

Introducing the artificial perturbation parameter £ and 
rewriting Eq.(4) yields the following iteration equation. 
£(u c ) tt (x,t) = 

(—2 + x + t(l + t(6 + 2f + t 3 + (2 + f 2 )x))j sinGc) 

= (1 + t 2 ) 3 

— (u n ) tt (x,t) + (t + x)(u n ) xx (x, t),n = 0,1,2,... (7) 
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Appropriate to the initial conditions, an initial estimation 
to begin the iteration process is proposed as u 0 (x, 1) = 
sinGc). Subrogating this initial condition in Eq.(4) and 
solving it gives 


(u c (x, 1))o : 


6-6 t-7V 


3 t 2 x — 3 t 4 x 


6(1 + 1 2 ) 

6(1 + t 2 )(l + tx)tan _1 (t) 
6(1 +1 2 ) 

3(1 + t 2 )(t — x)ln(l + t 2 ) 

6(1 + 1 2 ) 


sin(x) 

sin(x) 

sin(x) 


+ c t (x) + tc 2 (x ) (8) 

So the first iteration result using the initial conditions is 
u t (x, t) = u 0 (x, t) + e((u c (x, t)) 0 ) (9) 

or 


u^x, t) = sin(x) + 


6t 2 + (t + t 3 )( 6 +1 2 + 3 tx) 


6(1 + t 2 ) 

(1 + t 2 ) (2(1 + tx)tan _1 (t)) 

6(1 +t 2 ) 

(t- x) ln(l + t 2 ) . 

-sin ( x) 


sin(x) 
sin (x) 


6(.1+P) - SmW (10) 

If the procedure continues similarly, we get the following 

results 

u 2 (x,t) 

61 2 + (t + t 3 )(6 +t 2 + 3 tx) 


= — sin(x) + - 


6(1+ t 2 ) 


■ sin(x) 


3(1 + t 2 )( 2(1 + tx)tan 1 (t) + (t — x) ln(l + t 2 )) 


6(1 + 1 2 ) 


sin(x) 


r( 4C1 


720(1 + t 2 ) ' 

+ t 2 )cos(x) (t(—30 — 20t 2 + 9t 4 — 150tx + 15t 3 x) 

+ 30(t+ t 3 — x + 3t 2 x) ln(l + t 2 )) 

- (-1440 + 25201 3 + 4t 8 - 180t(-12 + x) + 24t 7 x 
+ 12t 5 (30 + 17x) - 30t 2 (—3 - 36x + 10x 2 ) 

+ t 6 (99 + 30x 2 ) — 5t 4 (—37 — 216x + 54x 2 ) 

- 30(1 + 1 2 )(-l + t 4 + t(36 - 8x) - 36x + 2x 2 

- 6t 2 x 2 ) ln(l + t 2 ))sin(x) 

+ 60(1 + t 2 )tan _1 (t)(—2(— 1 + t 4 — 6 tx + 2t 3 x) cos(x) 
+ (36 — 3x + 6 t 2 x + t 4 x + t(2 + 36x— 6x 2 ) 


+ 2t 3 (l + x 2 ))sin(x))) 


( 11 ) 


U3(X ’ 0 = 3175200(1+ t 2 ) X - 21(t + t 3 )C8400(-6+ x) 

+ t(t(—33600 + t(—15350 + 31(5040 
+ 608t+ 105t 3 ))) + 20t(—1778 
+ 31(420 — 3991 + 22t 3 ))x +630(54 

- 77t 2 + 2t 4 )x 2 - 180(27 

+ 1400x)))cos (x) + 3175200sin(x) + (1 
+ t 2 )(1260cos(x)(2(10t 7 + 70(—6 + x) 

- 210t 4 (—2 + x) - 420t 2 x + 56t 6 x 

+ 7t(—7 + 45(—8 + x)x) + 63t 5 (l + x 2 ) 

- 70t 3 (2 + 3x(—4 + 3x)))tan _1 (t) 

- (—17 — 840x + 7(2t 6 + 61(20 — 7x) 

+ 30t 5 x- 90t 2 (—4 + x)x + 9x 2 

+ 20t 3 (6 + x) + 5t 4 (5 + 9x 2 )))ln(l 
+ t 2 )) + (t(245t 8 + 2205t 7 x 
+ 1260(—47 + 42(—30 + x)x) 

+ 45t 6 (—155 + 126x 2 ) + 63t 5 (560 
+ 313x+ 70x 3 ) - 63t 4 (—9262 
+ 21x(—160 + 37x)) — 210t 2 (326 
+ 21x(—360+ 67x)) +6301(1260 
+ x(—1304 + 21x(—200 + 9x))) 

- 105t 3 (—7980 + x(—8966+ 21x(-120 
+ 77x)))) - 1260(—47 + 5t 7 x + 42(-30 
+ x)x — 315t 2 (—8 + x)x + 35t 4 (19 

+ 12x) + 21t 5 x(-2 + x 2 ) + 7t 6 (—10 
+ 3x 2 ) + 7t(120 +x(—104+ 15(—24 
+ x)x)) + 35t 3 (24 + x(35 - 6(-4 
+ x)x)))tan _1 — 630(—420 + 5t 7 + 71 6 x 
+ x(152 — 21(—40 + x)x) + 1051 2 x(—13 
+ 2(-12 + x)x) - 63l 5 (—6 + x 2 ) 

- 35l s (13 + 6x 2 ) + 35l 4 (12 + 14x 

- 3x 3 ) + 7l(—20 + 3x(—160+ 9x)))ln(l 

+ l 2 ))sin(x)) (12) 


hi this study, the Eq.(5) in [16] is solved by using PIA 
and the following first-order difference scheme 

z- 2 (u k+1 -2u k +u k ^) + A k u k+1 = f k , 

A k = MtfcX/fc = f{t k ),t k = kz, 

1 < k < N - 1, Nz = T, 

^t -1 (u 1 — Uq ) + i A^ 2 u t = iA^ 2 u 0 + xp,u 0 = cp. 
studied in [17] and the second-order difference scheme 

( U k +1 2u k +Ufc —1 + ll A 2 1J -f 

t 2 + 4 ^t“k+l —Ik’ 

. A k = A(t k ),fk = f(t k ),t k = kz,l <k<N- 1 ,Nt = T, stu 

(/ +T 2 A 0 )l -1 (u 1 — Uq) 

= \(fo-A qUo)+ xp,f 0 = f{0),u 0 = cp. 
died in [16]. The results are compared and discussed. 
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IV. FIGURES AND TABLES 



Fig. 1: Surface plot of the third order PIA solution. 


Fig. 2: Surface plot of the exact solution. 


Table. 1: Comparison of the third order PIA absolute 
errors 



Third order PIA Absolute Errors 

X 

t = 0.2 

t = 0.4 

t = 0.6 

27T/10 

3.26329E — 7 

2.32795 E 

- 5 

2.85770E 

- 4 

47T/10 

3.23876E — 7 

2.47931E 

- 5 

3.24057E 

- 4 

67T/10 

1.78211E — 64 

1.23072E 

- 4 

1.48412E 

- 3 


Table. 2: Comparison of the first order difference scheme 
absolute errors 



First Order Difference Scheme Absolute Errors 

X 

t = 0.2 

t = 0.4 

t = 0.6 

2n/10 

2.66413E — 4 

4.85849E 

- 4 

6.71907E 

_ 4 

47r/10 

4.33106E- 4 

8.00497E 

- 4 

1.12590E 

- 3 

67r/10 

4.35134E- 4 

8.14527E 

- 4 

1.16251E 

- 3 


Table. 3: Comparison of the second order difference 
scheme absolute errors. 



Second Order Difference Scheme 

X 

t = 0.2 

t = 0.4 

t = 0.6 

27r/10 

6.50953E — 8 

1.19193E 

- 7 

4.00733E 

- 8 

47r/10 

1.11504E- 7 

2.51531E 

- 7 

1.43050E 

- 7 

6 tt /10 

1.71319E— 7 

5.05788E 

- 7 

7.22128E 

- 7 


As shown in Table 1 and 3, the second order difference 
scheme is approximately 10" 3 times better than the first 
order difference scheme. On the other hand, the results 
obtained by PIA are better than the first-order difference 
method but they are not as satisfactory as the results 
obtained second-order difference method. 

For more steps of PIA, various partial differential 
equation scan be studied and solved and the results are 
compared with each other as future problems 

V. CONCLUSION 

In this paper the approximate solution of a wave partial 
differential equation is obtained by previously developed 
efficient method, perturbation-iteration algorithm The 
method gives highly approximate solutions after a few 
iterations. The results are compared with the exact 
solution via absolute error and finite difference method. 
For this purpose first and second order difference 
schemes are applied. Also some surface plots and tables 
are presented to show the reliability of the method. This 
confirms that the method is ready to apply for wider class 
of partial differential equations. 
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